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ABSTRACT 

We present the mid-infrared spectrum, obtained with the Spitzer Infrared Spectrograph (IRS), of 
HD 98800, a quadruple star system located in the 10-Myr-old TW Hydrae association. It has a known 
mid-infrared excess that arises from a circumbinary disk around the B components of the system. The 
IRS spectrum confirms that the disk around HD 98800 B displays no excess emission below about 
5.5 /im, implying an optically thick disk wall at 5.9 AU and an inner, cleared-out region; however, 
some optically thin dust, consisting mainly of 3-/im-sized silicate dust grains, orbits the binary in 
a ring between 1.5 and 2 AU. The peculiar structure and apparent lack of gas in the HD 98800 B 
disk suggests that this system is likely already at the debris disks stage, with a tidally truncated 
circumbinary disk of larger dust particles and an inner, second-generation dust ring, possibly held up 
by the resonances of a planet. The unusually large infrared excess can be explained by gravitational 
perturbations of the Aa-t-Ab pair puffing up the outer dust ring and causing frequent collisions among 
the larger particles. 

Subject headings: circumstellar matter — binaries: close — stars: individual (HD 98800) — planetary 
systems: formation — infrared: stars 



1. INTRODUCTION 

During the last few years, particular interest has arisen 
for a certain stage in the evolution of circumstellar disks 
surrounding T Tauri stars: so-called transition disks, 
which are characterized by cleared-out or optically thin 
inner disk regions and a truncated optically thick outer 
disk. They are thought to be in a ph ase when the disk 
dissipates rapid l y fro m inside out fe.g.. lStrom et al.lll989t 
ISkrutskie et al.l Il990 ). on a timescale of ~ 10^ years. 
A few disks in transition have been ide ntified in the 
1-2 Myr old Taurus star-forming region (iForrest et all 
[200l iD'Alessio et al.l[2005l iCalvet et al.ll2005D : the dif- 
ferent amount of material left over in the inner and outer 
disks of these transition objects implies different forma- 
tion m echanisms or stages in disk evolution (Calvet et al.l 
|2005| ). Planet formation, which is thought to occur in cir- 
cumstellar disks around T Tau ri stars, mi ght play a role 
in the cleari ng of i nner disks ((Marsh fc M ahonc v 119921: 
Calvet et al.1 120021: iQuillen et all 120041: ID'Alessio et all 
2005f ). but it is not well-known whether other mecha- 
nisms, such as ph otoevaporation, can produce these tr an- 
sition disks (e.g.. [Alexander. Clarke, fc Pringlell2006D . 

On the other hand, disks with inner gaps can also be 
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found in stable, long-lived configurations: a binary star 
system surrounded by a circumbinary disk might cre- 
ate an inner gap due to reso nant and tidal interactions 
(|Artvmowicz fc Lubowlll99"^ . thus generating a spectral 
energy distribution (SED) typical of transition disks. For 
example, the T Tauri star St 34, which is a spectroscopic 
binary, has an inner disk depleted in dust, likely an ef- 
fect of the gravitation al perturbations of th e binary on 
the inner disk regions (iHartmann et al.ll2005f l. Since this 
system is likely older than 10 Myr, its disk configuration 
cannot be attributed to a short-lived transitional stage. 
In the ~ 25-member TW Hya as- 

sociation, which is 5-15 Myr old 

(IStauffer. Hartmann. fc Barrado v Navascuei Il995t 
IWebb et all 119991 : IWeintraub. Kastner. fc Barvl l2000( l. 
four objects are characterized by significant infrared 
excesses: T W Hya, Hen 3- 600, HD 98800, and HR 
4796A (e.g.. iLow et al.ll2005h . Except for HR 4796A, 
whose weaker IR excess clearly places it into the debris 
disk category, these stars are surrounded by substantial, 
probably protoplanetary, disks. Both TW Hya and Hen 
3-600 display close to no disk emission below 7 fim, 
an in dication of truncated inner disks (jUchida et al.l 
l2004f ). In particular, TW Hya has been suggested as 
the forma tion site of a prot oplanet based on its inner 
disk gap (ICalvet et al.l 1200^ . Hen 3-600 consists of a 
spectroscopic binary (the A component) and a com- 
panion separated by 1'.'4 (e.g., IWebb et al. 199 9;); the 
disk s urrounds only Hen 3-600 A (jJavawardhana et al.l 
Il999a[ ). implying that the inner disk gap might be a 
result of gravitational perturbations as in the case of St 
34. 

In this paper we introduce the mid-infrared spectrum 
from 5 to 36 fiui of HD 98800 ( TV Crt) obtained w ith 
the Infrared Spectrograph*^ (IRS: iHouck et al.ll200^ on 

* The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
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board the Spitzer Space Telescope (jWerner et al.|[200^ . 
HD 98800 is a quadruple system consisting of a visual 
binary with a projected separation of 0'.'8 (which corre- 
sponds to 38 AU at the distance of 47 pc determined 
by Hipparcos) whose components are spe ctroscopic bina - 
ries with separations of about 1 AU (fBoden et al.l l2005l) . 
Mid-infrared imaging revealed that the strong infrared 
excess of the system arises from component B only; in ad- 
dition, it sets in at about 7 iim, implying cleared-out in- 
ner disk regions (jGehrz et al.lll999l: iKoerner et al. ''2000: 
iPrato et"alll2001h . Our mid-infrared spectrum confirms 
that the disk around HD 98800 is a transition disk and 
allows us to determine the location of the circumstellar 
dust, as well as to derive its mineralogical composition. 
This paper is structured as follows: in § 2 we present our 
observations and data reduction; in § 3 we construct the 
SED of HD 98800 B, and fit dust and disk models to it; 
and in § 4 we present a discussion of this transition disk 
and our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

HD 98800 was observed during IRS campaigns 9 and 
17 on 2004 June 25 and on 2005 January 03, respec- 
tively. The observations of campaign 17 repeated those 
of campaign 9, since the observations of the earlier cam- 
paign were somewhat compromised by a bright satura- 
tion event. However, the observations of campaign 17 
were slightly mispointed; since we were able to miti- 
gate the effects of the saturation on the array (apparent 
mostly in the form of additional rogue pixels), we used 
the IRS spectrum of campaign 9 for our analysis. 

To obtain the full mid-IR spectrum, we used the Short- 
Low (SL; 5.2-14 /im; A/AA=60-120), Short-High (SH; 
9.9-19.6 Aim; A/AA==600), and Long-High (LH; 18.7-37.2 
/im; A/AA=600) modules. For a clearer presentation of 
the IRS spectrum, we rebinned the SH and LH spectra 
to a resolution of 300, and we truncated the SH spec- 
trum below 14 /im. The observations were carried out in 
mapping mode, where we mapped the target in 3 steps 
separated by three-quarters (for SL) or half (for SH and 
LH) of the slit width in the dispersion direction and 2 
steps separated by a third of a slit length in the spatial 
direction. 

We extracte d and calibrated our data with the SMART 
software tool ijHigdon et al.|[200^ after fixing bad pix- 
els in the arrays using a simple interpolation over good, 
neighboring pixels in the spectral direction. In SH and 
LH, we also fixed all so-called "rogue" pixels identified 
in darks from campaigns 1 to 18 (for LH) or 1 to 28 (for 
SH). We subtracted the background of our SL spectrum 
by using the observation taken in the other nod position 
in the same order, then extracted the spectrum using 
a variable-width column extraction that scales with the 
width of the point-spread function. We were not able 
to subtract any background in our SH and LH spectra, 
which we extracted with a full-slit extraction. The SL 
spectrum was calibrated with a Lac (Al V), the SH and 
LH spect ra with £, Dra (K2 III), using template spec- 
tra from iCohen et al.l (|2003[ ) . The final spectrum was 
obtained by averaging the spectra of the two central 
map positions. Given that the SL, SH, and LH spec- 
tra stitched together without applying any scale factors. 

Jet Propulsion Laboratory and the Ames Research Center. 
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Fig. 1.— SED of HD 98800, composed of the IRS spectrum 
(which includes t he flux of all four stars in the system), B, V, I 
photometry from lSoderblom et al.1 i ITQQI). the 2MASS J , H, and 
Ks fluxes, M, TV, Q photometry from IPrato et al.l l|2001l '). and the 
IRAS 12, 25, 60, and 100 fj.m fluxes. The photosphere of a K5 
star ( based on the photospheric colors from lKenvon fc HartmannI 
l|1995l ) and normalized at J) is also sketched in. All fluxes were 
dereddened using Av=0.44 and Mathis's reddening law (MathlJ 
[19901) . 

the background contributions in SH and LH are likely 
negligible. 

By comparing the 10 and 24 /xm emission derived 
from our IRS spectrum of HD 98800 with ground- and 
space-based measurements (iJava wardhana et al.lll999bl : 
[Low et al. 2005), we infer an absolute spectrophotomet- 
ric accuracy of ~ 5%. Our relative accuracy, determined 
by the scattering of neighboring fiux values, is higher in 
SL than SH and LH: features above the noise level are 
very likely real in SL, while LH is still dominated by 
calibration artifacts. 

3. ANALYSIS 

3.1. Spectral Energy Distribution 

Previous observations of HD 98800 indicated that 
the infrared excess originates from a circumbinary disk 
around the component B spectroscopic binary. The SED 
is consistent with it being a trans ition disk, extending 
between a few an d ~ 10-20 AU (jKoerner et al.l l2000l : 
IPrato et al.ll200il ). Observations of the Ha line, which 
did not distinguish the emission from the four compo- 
nents, revealed a total Ha equivalent width of 0.9 A and 
a width comparable to that of Pleiades K dwarfs, sug- 
gesting an origi n in chromos pheric activity and not in an 
accretion fiow (jSoderblom et al.lll996l ). An upper limit 
to the mass accretion rate has not been determined, but 
it is likely well below 10~^^ AIq yr~^ (|Muzerolle et al.l 
I2000f) : thus, accretion onto HD 98800 B has virtually 
come to a halt. The absence of accreti on signatures 
(jSoderblom et all [199^ iWebb et all Il999h is in accor- 
dance with the lack of dust emission from the inner disk, 
i.e., the outer disk is prevented from accreting toward 
the star. 

Our Spitzer IRS observations included emission of all 
four components of the system, since the narrowest IRS 
slit is 3'.'6 wide. Therefore, in the SED shown in Fig- 
ure [H we used unresolved optical to far-IR photome- 
try of HD 98800 (see figure caption for details) and the 
IRS s pectrum, all corrected for reddening using ^^=0.44 
(e.g., jS oderblo m et al.l ll998) and Mathis's reddening law 
(jMathis J^90. )- Also shown is the photospheric emissio n 
of a K5 star (jSoderblom et al.lll998t IPrato et al.ll2001h . 
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Fig. 2.— SED of HD 98800 B, composed of optical to mid-IR pho- 
tometry compiled by Prato ct al. (2 00l|) (diamonds), IRAS fluxes 
(squares), and the MIPS fluxes of I Low et al.l 120051) (triangles). 
The photosphere is also sketched in, based on the photospheric 
colors of a K5 star IjKenvon &: Hartmann' 119951) and normalized 
at the J-band flux given in [Prato ct al. (2001^). The IRS spectrum 
was corrected for the contribution of the A component (see text for 
details). All fl uxes wer e dcrcddened using ^4^=0. 44 and Mathis's 
reddening law l lMathislll99Q') . 

based on the photospheric co lors of a star with this spec- 
tral type (iKenvon &: Hartma nn 1995) and normalized at 
the unresolved 2MASS J-band flux. Despite the fact that 
our IRS observation does not distinguish the emission 
from A and B, the SED suggests that none of the four 
components generates excess emission above the level ex- 
pected from a stellar photosphere out to about 6 /im. 
This supports the idea that the inner disk of HD 98800 
has been cleared out. 

In order to retrieve the true IRS spectrum of the B 
component, we estimated the contribution of the A com- 
ponent and subtracted it from the spectrum. Com- 
ponent A does not dis play an infrar ed exce ss and has 
an extinction Ay ~ ([Koerner et al. 2000; Prato et al.l 
1200 It iBoden et all 120051), while Ba-hBb seem to lie be- 
hind Ay =0.44 (e.g., [Sodirbiomiianillil). This dif- 
ference in extinction likely indicates that we are ob- 
serving HD 98800 B through some circumbinary dust, 
as first noted by Jbkovinin (1999); this could also 
explain the slight pho tometric variability observed by 
ISoderblom et al.l (|1998D . 

We estimated the mid-IR emission of HD 98800 
A by using the p hotospheric colors of a K5 star 
(jKenvon fc Hartman n 1995) norma lized at the J-ban d 
flux of the A component given in iPrato et al.l (|2001h . 
This estima t e is c onsistent with the measurements by 
iPrato et"al] (|2001h . who determined that the contribu- 
tion of A to the system flux decreases with wavelength 
(50%, 44%, 30%, and 20% at 2.2, 4.8, 7.9 ^i^Jl, and 
8.8 /im, respectively, and less than 10% at longer wave- 
lengths). The subtraction of the estimated emission of 
the A component from the IRS spectrum will affect both 
the flux level and shape of the spectrum particularly in 
the 5-8 fim wavelength range; given that we do not ac- 
count for atmospheric features, the resulting spectrum 
carries a higher uncertainty in this wavelength region. 

The SED of HD 98800 B is shown in Figure El it was 
constru cted with photome try of the B componen t com- 
piled b y IPrato et all (j2001h . MIPS fluxes from Lo w et all 
^2005^, IRAS Ruxes, the IRS spectr um corrected for the 
contribution of A, and the photosphere of B (approxi- 
mated using the colors of a K5 star normalized at the 



J-band flux of B) . The emission of B is photospheric out 
to about 5.5 /xm; beyond that, the infrared excess sets 
in rather sharply. The steep rise of the SED at 8 /im is 
reminiscent of that observed in CoKu Tau/4, a transi- 
tion disk in the Taur us star-forming r egion with an inner 
disk hole of 10 AU (jPorrest et all [2004 : D'Alessio et all 
l2005h . The SED peaks at 20 /zm, then decreases as a 
single-temperature blackbody. The outer disk is likely 
truncated by component A, who se distance from B is al- 
most 50 AU at closest approach (jTorres et al.l ll995!: note 
that the orbit of B around A is very eccentric and seen al- 
most edge-on). Thus, the disk around the spectroscopic 
binary HD 98800 B seems to be truncated both at the 
inner and at the outer disk edge. 

3.2. Dust Composition 

iKoerner et al.l ()2000l ) inferred from narrow-band pho- 
tometry in the 10-/im region that the silicate emission 
feature of HD 98800 B is broader and more structured 
than expected from amorphous silicates, indicating the 
presence of crystalline silicates. The IRS spectrum al- 
lows us to analyze the 10-/im silicate feature in greater 
detail and to derive the dust species that generate the 
infrared excess emission. 

First, we subtracted the photospheric emission, con- 
structed from the colors of a K5 star normalized at the 
J-band flux, from the dereddened IRS spectrum of HD 
98800 B. Next, given the large infrared excess and the 
presence of 10 and 18 /im silicate emission features, we 
fit the residuals with two components, an optically thick 
region modeled as a blackbody, and a warmer, optically 
thin region of constant source function (also modeled as 
a blackbody): 



e-^-), (1) 



where Tthick is the blackbody temperature of the opti- 
cally thick region, Tthin is that of the optically thin re- 
gion, is the wavelength-dependent optical depth of the 
optically thin cloud, and Qthick and Qthin are the solid 
angles of the emitting regions for the optically thick and 
thin components, respectively. 

The optically thick component can be thought of as 
consisting of large (mm-sized) dust grains, which do not 
generate spectral features and whose opacities are inde- 
pendent of wavelength, thus creating a continuum com- 
ponent with a blackbody-type emission. Alternatively, a 
very large mass of small grains will generate blackbody- 
shaped, optically thick emission. Thus, the detailed com- 
position of the dust only enters in the optically thin com- 
ponent. 

For optically thin emission, the fiux becomes linearly 
dependent on the optical depth, which can be expressed 
as the sum of the optical depth of each dust component 
= K^iNi, where k^- and Ni are the mass absorption 
coefficient and column density for each dust species, re- 
spectively. Since the mass of each dust species is derived 
via rrii — NiQ,thin<P' , where within is the solid angle sub- 
tended by the dust and d is the distance to the system 
(assumed to be 47 pc), the total optical depth of the 
optically thin component can be expressed as 
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TABLE 1 

Dust Component Fits of HD 98800 B 



Model 


Dust Species 




thin 










[K] 


m 


[10-" sr] 
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5 fim amorphous pyroxene 


155 


373 


1.39 


18.6 
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small amorphous carbon 












+ 5 iim amorphous pyroxene 


151 
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1.47 


12.1 
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small amorphous carbon 












+ 3 fim amorphous pyroxene 


153 


298 
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Fig. 3. — Dust model fit, using ISM-type dust grains, for HD 
98800 B {thick gray line), whose IRS spectrum is plotted after 
photosphere subtraction; note that the error bars for the IRS data 
points are generally smaller than the plotting symbols and are thus 
not shown. The model components are a blackbody at temperature 
Tss=161 K for the optically thick part {dash-dotted line) and 
sub-/im amorphous olivine at T=662 K for the optically thin part 
{dashed line). The solid angle of the emitting area of the optically 
thick component is equal to 1.23 X 10" sr. The large value 
of 55.4 (calculated between 6.5 fira and 1.3 mm) proves that this 
model is not a good fit. 

where the sum is carried out over the number N of dust 
components. Thus, the dust emission is calculated as 

1 ^ 
Fy = ^thickBi,{Tthick) + -^By{Tthin) ^ Huirrii (3) 

1=1 

We solved the linear model for the best-fitting set of tem- 
peratures Tthick and Tthm-, solid angle Q.thick-, and dust 
masses mi, starting with just one dust component and 
adding components as necessary to improve the fit. 

For the optically thin dust species, we adopted 
amorphous carbon and silicates, which are typical for 
the i nterstellar me d ium. We used optical constants 
from [ Draine k Led (Il98l for A < 7.5 /im and from 
iDorschner et al.l (llQQSlT for A > 7.5 /im (using power 
laws to extrapolate values from 200 /im to 2 mm) for 
amorphous silicates of olivine (MgFeSi04) and pyrox- 
ene ( Mgo.sFep. 28103) com position, and optical constants 
from IZubko et al.l (|1996) for amorphous carbon. For 
sub-micron grains, we used a CDE2 shape distribution 
([Fabian et al.l l200ll ) to compute opacities, while larger 
grains (> 1 /xm) were adopted to be porous spheres with 
a 50% volume fraction of vacuum. 

The first model attempts, using only sub-micron (i.e., 
ISM-like) olivine or pyroxene grains as the optically thin 



component, yielded poor fits (for an example, see Figure 
[3]) ; the emission at the long- wavelength side of both the 
10 and 18 /.tm silicate feature is underestimated, suggest- 
ing the presence of larger grains. Amorphous pyroxene 
represents a better fit to the observed peak positions of 
the 10 and 18 /xm features than amorphous olivine, so it 
was assumed as the dominant large grain species. 

The best-fitting models, ordered by decreasing reduced 
values (x^), which were calculated from 6.5 ^m to 1.3 
mm, are shown in Figure S] with the residuals of HD 
98800 B. The adopted model components and their pa- 
rameters are displayed in Table [T] We note that most of 
the structure seen in the 5-6.5 /im region of the residuals 
is likely due to photospheric features from the four stellar 
components of the system, which we did not account for 
in our photosphere subtraction process. Thus, we expect 
a mismatch between models and residuals in this part of 
the spectrum, and we did not include these data points 
in the computation of our xj values. 

We first tried a model with an optically thick com- 
ponent at 155 K and a single optically thin component 
consisting of 5 /um amorphous pyroxene (Fig. |4]a). This 
fit improved by adding sub-micron amorphous carbon 
to the dust mixture (Fig. |3]b); in particular, the short- 
wavelength excess starting at about 6.5 /xm is fit better. 
A slightly better fit of the silicate feature was achieved 
by using 3 /im amorphous pyroxene and small amorphous 
carbon as the optically thin component (Fig. [5]c). Fi- 
nally, the best fit was obtained by adding 3 ^m. amor- 
phous olivine to the previous mixture and adopting an 
optically thick component at 152 K (Fig. |l]d). For all 
of the models shown in Figured! the adopted solid angle 
lies in the (1.4-1.5) x 10~^^ sr range, and the temper- 
ature of the optically thin grains is higher than that of 
the optically thick component by about a factor of two 
(see Table [iJ. Even though the differences between these 
four models are small (especially if judged by eye), the 
best fit determined here for the optically thin disk region 
around HD 98800 B will be confirmed by our disk model 
in §[S31 

The best model fit and its components are shown in 
Figure [5] with the residuals of HD 98800 B over the mid- 
infrared to mm wavelength range. The model fits the en- 
tire excess emission from 6.5 /im to 1.3 mm remarkably 
well. Table [2] lists the parameter values of this model 
together with their l-cr uncertainties, which were esti- 
mated from models lying within Ax^ = 1 of the best- 
fitting model. Since the optically thick component is 
more dominant and tightly constrained by all the mea- 
surements beyond about 14 /im, the relative uncertain- 
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Fig. 4. — Best dust model fits for HD 98800 B (thick gray lines); the IRS spectrum is plotted as in Figure |3] The values for each model 
are also shown in each panel. The fits consist of blackbodies at different temperatures (Tbb) for the optically thick component [dash- dotted 
lines) and different dust species for the optically thin component (dashed lines): (a) Tss=155 K, and 5-/im-sized amorphous pyroxene 
grains; (b) Tgg=l51 K, and small amorphous carbon and 5-/im-sized amorphous pyroxene grains; (c) Tgg=153 K, and small amorphous 
carbon and 3-/im-sized amorphous pyroxene grains; (d) TsB=f52 K and small amorphous carbon, 3-/im-sized amorphous pyroxene, and 
3-/im-sized amorphous olivine grains. 

Hydrae association, TW Hya (K7 spectral type) and Hen 
3-600 A (M3 spectral type), are dominated by smaller 
grains and are characterized, in particular in the case of 
Hen 3-600 A, by a n otable fraction of crystalline grains 
(jSargent et al.l |2006f). Our dust model fit to HD 98800 B 
excludes any sizable amount of small crystalline silicates, 
suggesting that the la rger width of the 10 li ra silicate fea- 
ture already noted bv lKoerner et all (|2000f) should likely 
be attributed to grain growth alone. 



TABLE 2 

Parameters and Uncertainties for the Best Dust 
Component Fit of HD 98800 B 



Parameter 



Value 



^pyroxene 
^olivine 
^carbon 



1.46 ±0.02 X 10"^^ sr 
151.8 ±0.4 K 
310±2 K 

2.3 ± 0.2 X lO""* lunar masses 
1.7 ±0.2 X lO"** lunar masses 
6.9 ± 0.8 X lO-'^ lunar masses 



Note. — The mass of the moon is equal to 7.35 x 10 g. 

ties for both ^tkick and Tthick are smaller than for the 
parameters describing the optically thin component. We 
derive a combined mass for the optically thin dust grains 
of 4.7 X 10^'' lunar masses. The mass of the optically 
thick component cannot be constrained, but it is likely 
several orders of magnitude larger than the mass of the 
smaller grains. 

The fact that dust grains larger than ISM-type amor- 
phous grains (whose sizes are < 1 /im) constitute the 
bulk of the optically thin component indicates that dust 
growth must have occurred in the disk of HD 98800 B. 
The other two objects with substantial disks in the TW 



3.3. Disk Model 

In order to derive the structure of the circumstellar 
material of HD 98800 B, in addition to its composi- 
tion, we compu t ed a mode l following the met hods of 
iD'Alessio etall (I20Q1 and iCalvet et al.l (l200l . The 
stellar luminosity of the Ba+Bb pair was adopted to be 
0.7 L0, a value obtained by integrating the photospheric 
fluxes of the unresolved B component, and the effective 
temp erature was assumed to be 4 350 K, typical for a K5 
star (|Kenvon &: HartmannI [19951 ) . The model has two 
components: an optically thick wall and an optically thin 
region, corresponding to the blackbody component and 
the optically thin dust grains, respectively, from the dust 
component fit of § 13.21 Even though the model implicitly 
assumes the presence of gas in addition to the dust, no 
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Fig. 5. — The best-fi tting dust model of HD 98800 B, plotted with t he photosph ere-subtracted IRS spectrum, long-wavelength 
photometry compiled by IPrato et al.l pOOll ). MIPS 70 and 160 ^tm fluxes l|Low et al.l 12005). and SCUBA 1350 and 2000 pira fluxes 
ilSylyester. Dunkin. fc Barlowll20'0lfr The various model components are listed in the top right quadrant of the figure. 



assumptions regarding the gas were made, since it does 
not enter any of the calculations. 

The wall emission was calculated assuming amorphous 
carbon and silicate grains with sizes between 1 and 3 
mm. Since the grains are so large, their opacities are 
gray (i.e., independent of wavelength), and thus their 
composition does not affect the fit; their emission cor- 
responds to that of a scaled blackbody. Therefore, this 
wall matches the optically thick component of § 13.21 The 
composition of the optically thin component was adopted 
from § 13.21 namely sub-micron amorphous carbon, 3 
^m-sized amorphous olivine and 3 /im-sized amorphous 
pyroxene grains. The small amorphous carbon grains, 
which absorb stellar radiation more efficiently than the 
larger grains, are necessary to heat the dust and yield 
enough emission in the 6-8 /im range; the 3 /im-sized 
grains generate emission that provides a good fit to the 
silicate features at 10 and 20 //m. 

The results of the model calculations are shown in Fig- 
ure [6] and in Table [3l We show different model fits which 
result from varying the contribution of the optically thin 
and thick components. As in the previous section, the 
properties of the inner disk wall are well-constrained by 
the flux measurements beyond about 14 /xm; in order to 
obtain the correct blackbody shape, it has to be placed 
at a distance of 5.9 AU from the B pair, where the tem- 
perature in its upper layers (closest to the star) amounts 
to 140 K. Not surprisingly, this temperature is very sim- 
ilar to the one derived in § 13.21 for the optically thick 
component. The height of the wall above the midplane 
is treated as a scale factor for the wall emission; it is 



adjusted such that the emission from the optically thick 
and thin components yields the best possible fit. For 
the best-fit model (Figure [6] d) , this height is 0.75 AU. 
The solid angle subtended by the wall amounts to 1.57 x 
2Q-13 gj. (assuming a distance of 47 pc and an inc lination 
angle of 67°, as determined bv lBoden et al.|[2005f ). which 
is just somewhat larger than the value of 1.46 x 10^^^ 
sr derived from the best-fitting dust model in § 13.21 

The optically thin component has to lie inside the re- 
gion delimited by the optically thick wall in order to pro- 
vide enough emission at wavelengths between 5 and 20 
/im; in addition, its innermost location and radial extent 
are fairly well constrained by the observed 10-20 /im flux 
ratio. As can be seen in Figure [6] a, an optically thin re- 
gion extending from 2 to 5 AU will underestimate both 
the short- and long-wavelength excess emission; decreas- 
ing the outer radius results in a better match of the flux 
beyond 20 /im (Fig. [6] b) . Concentrating the optically 
thin region close to 2 AU yields a very good fit (Fig. [B] 
c), while the best fit to the observations is obtained by 
an optically thin ring located between 1.5 and 2 AU (Fig. 
[6]d). The dust grains in this ring are at a temperature 
between 307 and 266 K, similar to the value found for 
the best-fitting dust model of § 13.21 

The wall radius of 5.9 AU we derive is comparable 
to the inner radius ran ge of 5.0 ± 2.5 AU determined 
bv lKoerner et al.) (|2000D from simple SED models; both 
values are l a rger t han the estimate of 2 AU derived by 
iPrato et al.l ()2001[ ) from fits to their unresolved mid- 
infrared ima ges. However, our result is consistent with 
the fact that IPrato et al.l (|2001h did not resolve the disk: 
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Fig. 6. — Sequence of disk models for HD 98800 B (thick gray lines); the SED is plotted as in Figure[2] Also shown are the components 
of the disk models, the optically thick disk wall (dash-dotted lines) and the optically thin inner disk region (long-dashed lines). The four 
panels display different models where the optically thin region consists of a ring between ri„ and rout, and the optically thick wall at 5.9 
AU is scaled by means of the wall height to result in a good overall fit (see Table [3] for details) . 



TABLE 3 
Disk Model Fits of HD 98800 B 



Model Optically thin component Optically thick component 

rin rout no T(rin) T(rout) Rwall H^aii 

[AU] [AU] [K] [K] [AU] [AU] 



a 2 5 0.04 266 175 5.9 0.56 

b 2 4 0.04 266 193 5.9 0.68 

c 2 2.2 0.22 267 255 5.9 0.75 

d 1.5 2 0.06 307 266 5.9 0.75 



Note. — no is used a scaling factor for the emission of the optically thin component: Fv = B^(Tdust)Ku-^^Y^^^2'Krdr 



the optically thin region is very compact and close to the 
binary, and even though at 12 /xm the brightness of this 
region and that of the wall are comparable, the large 
inclination of the wall might have prevented its direct 
detection. 

The optical depth of the optically thin component is 
displayed in Figure O it was computed using the best- 
fitting dust component model from § 13.21 and adopting 
a ring between 1.5 and 2 AU to determine ^thim which 
amounts to 5.85 xl0~^^ sr. Shown in the figure are the 
vertical optical depth [solid line) and the optical depth 
along the line of sight (^dashed line), assuming an inclina- 
tion angle of 67°. The latter quantity is actually an upper 
limit to the optical depth along our line of sight, since 
the thickness of the optically thin region is not known. 



Thus, since r^, lies between 0.02 and at most 0.13 over 
the mid-infrared range and decreases steadily towards 
longer wavelengths, the inner disk region will appear as 
optically thin from infrared to millimeter wavelengths. 

A sketch of the HD 98800 system is shown in Figure [8l 
There are two rings of material around the B component: 
an optically thick wall at 5.9 AU with a small radial 
extent, and an optically thin inner region between 1.5 
and 2 AU, just outside the binary orbit. We will discuss 
the implications of this peculiar structure in the next 
section. 

4. DISCUSSION AND CONCLUSIONS 

The mid-infrared spectrum of HD 98800 B reveals that 
it is surrounded by a transition disk whose infrared ex- 
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Fig. 7. — Optical depth as a function of wavelength for the op- 
tically thin component of the best-fitting dust model of HD 98800 
B. The solid line shows the vertical optical depth, while the dashed 
line represents an upper limit to the optical depth along the line 
of sight. 

HD 98800 A 



HD 98800 B 

wall 




optically thin dust 

Fig. 8. — Sketch of the HD 98800 system. Note: Distances and 
sizes of the disk components around HD 98800 B are not plotted 
to scale. 

cess emission starts at about 5.5 /im and decreases be- 
yond 20 /im. Its structure is some what reminiscent o f 
that of the tr ansition disks TW H ya (|Calvet et al.ll2002f l 
and GM Aur (jCalvet et al.ll2005f) . which also harbor op- 
tically thin regions inside an outer, optically thick disk. 
However, as opposed to the these two objects, HD 98800 
is a multiple system: the Ba and Bb stars form a close 
(~ 1 AU) eccentric (e = 0.78: lBoden et al.ll2005D binary 
system, and in addition the A pair comes as close as 50 
AU to the B pair, with the orbital planes of the two bi- 
naries inclined by a bout 30° wit h resp ect to each other 
([Torres et afl 119951: [Boden et alll2005f) . Thus, gravita- 
tional perturbations likely play a role in the confinement 
and distribution of the dust around HD 98800 B. 
A circumbinary disk around an eccentric binary is ex- 

S iected to be tidally trunca ted at the inner disk edge 
Artvmowicz fc Lubowlll99^ . If the eccentricity is large, 
several higher-order resonances play a role in truncat- 
ing the circumbinary disk by exerting tidal torques fur- 
ther out in the disk; thus, gap sizes of ~ 3.5a (where 
g is the semimajor axis of th e binary) can be explained 
(| Artvmowicz fc Lubowl [T994f ). However, these gravita- 
tional interactions apply to a gaseous disk, and the 
amount of gas in HD 98800 appears to be very small: 
the 1- cr upper limit placed bv lDent. Greaves, fc CoulsonI 
()2005f l on the sub-mm ^^CO J=3-2 line indicates an up- 
per limit of ~ 4 X 10"'' Jupiter masses (which is less 
than 1/24 of the gas mass of TW Hya), assuming opti- 
cally thin ^^CO emission. Some warm gas could still be 



present in the inner disk regions, since the sub-mm ^^CO 
line traces cold molecular gas in outer disk regions (where 
temperatures are bel ow ~ 100 K), but give i i the lack of 
accre tion signatures (jSoderblom et al.lll996l : I Webb et al.l 
Il999| ). the amount of gas in the inner disk should be 
small. 

In the absence of gas, the truncation of the outer 
disk at 5.9 AU could still be explained by the effect 
of resonances. Using the em pirical formula determined 
bv iHohnan fc WiegertI (|1999), who studied stable parti- 
cle orbits around an eccentric binary, the smallest, sta- 
ble orbit around HD 98800 B (which has a semimajor 
axis a ~ 0.98 AU, eccentricity e = 0.78, and mass ra- 
tio n = M2/(Mi + Ma) = 0.45; iBoden et al.ll2005[ ) is 
at 4.06 AU, somewhat smaller than the result we ob- 
tained by fitting a disk model to the SED. However, 
iHolman fc WiegerB ()1999[ ) note that the transition be- 
tween stable and unstable orbits is likely not sharp due 
to the effects of overlapping mean motion resonances. 
They also found indications that, over time, the stable 
region moves further out for high binary eccentricities. 

While a tidally truncated disk at > 4 AU could ex- 
plain the location of the disk wall, it would not be able 
to account for the presence of optically thin dust grains 
at 1.5-2 AU. The most likely interpretation for the un- 
usual structure of the circumstellar material of HD 98800 
B, also considering the probable scarcity of gas, is that 
the disk is already at the debris disks stage, when dust 
is its main constituent and is generated by collisions of 
larger bodies. The dust would be second-generation dust 
and not primordial material that survived for 10 Myr. 
The sub-micron carbon and 3-/xm-sized silicate grains in 
the disk would be replenished in planetesimal collisions, 
which might occur in the optically thick ring at ~ 6 AU. 
Due to Poynting-Robertson (PR) drag, the larger dust 
grains would spiral in from 6 AU towards the binary 
on a timesca, l e of a few 10^ years; since tpR <x. aD'^ /L,, 
(jBurns et al.lll979f) . where a is the grain size, D = 6 
AU and = 0.7 Lq, smaller grains would migrate even 
faster. Since HD 98800 is about 10 Myr old, the opti- 
cally thin, inner ring must be continuously replenished, 
assuming it is a long-lived structure. 

If the inner dust ring is explained by PR drag, then the 
absence of grains from 2 to 5.9 AU is puzzhng; drag forces 
acting on dust grains should distribute the dust uni- 
formly inside the radius at which the dust grains are cre- 
ated. A possible explanation for the observed gap could 
be a planet that formed just outside the unstable region, 
i.e., close to the inner disk wall, and that is temporarily 



at one of its inner mean motion resonances 1 


Liou fc Zoold 


ll997HMoro-Martm. Wolf, fc Malhortall2005 


). In the case 



of TW Hya, a planet was t hought to be respo nsible for 
clearing out the inner disk (jCalvet et al.ll2002f ). so it is 
conceivable that a planet also formed in the roughly co- 
eval HD 98800 system. On the other hand, this system 
is probably governed by complex dynamics due to the 
presence of four stellar components, implying overlap- 
ping resonances and variable gravitational perturbations. 

Given that the evidence supports that the disk around 
HD 98800 B is rather a debris than a protoplanetary 
disk, the presence of an optically thick dust component 
and the large infrared excess (L/fl/LboZ =17%) of HD 
98800 B seem unusual, but they might be explained 
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by the gravitational perturbations of the Aa+Ab pair. 
This type of perturbation can pump up eccentricities 
and inchnations of particles, and cause particles to be 
trapped in mean motion resonances, as was likely the 
case for Kuiper Belt objects under the influence of the 
giant planets and possibly a clo se encounter by a pass- 
ing, nearby star (e.g., Duncan, Levison. fc Buddl 119951 : 
llda. Larwood. fc Burkert. ,2000: iGladmanI I2005D . Peri- 
odic stirring of planetesimals in the outer disk around 
HD 98800 B by the A pair could be responsible for gen- 
eratin g copious amounts of dust (see iKenvon fc BromlevI 
|2002[ ). HD 98800 B is thus a unique type of debris disk, 
whose infrared excess is elevated to levels comparable to 
that of protoplanetary disks due to the particular config- 
uration of the four components in this system, resulting 
in gravitational perturbations that prevent the dust from 
settling into a flat disk. 

Even though HD 98800 appears to be a very dynam- 
ical system, it is unlikely to be in a short-lived transi- 
tional stage with ongoing clearing processes, in which 
the outer disk is being progressively eroded. HD 98800 
B belongs to a similar class of transition disks as St 34 
and Hen 3-600 A: a tight binary is responsible for tidal 
and resonant interactions with the disk, thus creating 
a stable, tidally truncated circumbinary disk. In addi- 
tion, collisions between planetesimals in the outer disk 
of HD 98800 B cause a collisional cascade of smaller 
grains, which then migrate towards the central binary. 
This outer disk is likely truncated due to the presence 
of the other component, as is the case with Hen 3-600 
A (fjavawardhana ct al.| Il999a[ ). However, St 34 and 
Hen 3-600 A are still accreting material, albeit at low 
levels (2.5 x 10"^° Mq yr y^ for St 34, ~ 5 x 10"" 
Ma, yr-^ for Hen 3- 600 A; iWhite fc Hillenbrandl 120051 : 
iMuzerolle et al.ll2000f) . while HD 98800 B is probably not 
accreting any more. Thus, St 34 and Hen 3-600 A are 
likely surrounded by evolved protoplanetary disks, while 
the disk around HD 98800 seems to have evolved even 
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further. The similarities and differences between the HD 
98800 and Hen 3-600 systems could increase our under- 
standing of disk evolution; the larger separation between 
Hen 3-600 Aa-|-Ab and B and the fact that all three stars 
are of later spectral type might play a role in the longer 
survival of primordial disk material around this system. 

The analysis presented in this paper suggests that 
transition disks comprise a varied group of objects, 
where different processes are responsible for creating a 
"transition disk" appearance, and where the timescales 
involved can vary substantially. Increasing the sample of 
transition disks from various star-forming environments 
and at different ages will shed light on the processes 
sculpting the disks. 
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